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1. Introduction 



Understanding the large distance behavior of asymptotically free gauge theories in terms of 
the weakly coupled short distance degrees of freedom is a major challenge for theoretical physics. 
In pure gauge theory with the standard Wilson's action, the available numerical data on L!^ lattices 
indicates that there is no phase transition for SU (2) or SU (3) and the theory should be in the 
confining phase for all values of the coupling. Convincing arguments have been given ^ in favor 
of the smoothness of the renormalization group flows between the two fixed points corresponding 
to the two limits. This suggests that it is possible to match the weak coupling and the strong 
coupling expansions of the lattice formulation. However, if we consider the two expansions, for 
instance for the average SU{2) plaquette as a function of j8 = there is a crossover region 
(approximately 1.5 < j8 < 2.5) where none of the two expansions seem to work. This situation 
can probably be explained in terms of singularities in the complex j8 plane that at this point 
are not completely understood. In these proceedings, we discuss the weak and strong coupling 
expansions of the density of states for SU (2) and compare them to Monte Carlo calculations. The 
density of states is the inverse Laplace (or Borel) transform of the partition function. Its logarithm 
can be interpreted as a "color entropy". This is discussed in section ^ where the basic concepts are 
defined. 

For the one plaquette model, the density of states is a function that has better convergence 
properties than the partition function [||]. This is explained in section ^. We would like to know 
if this property persists on lattices. The comparison between weak and strong expansions and 
numerical calculations of the density of states for a 6"* lattice are summarized in section ^ More 
details can be found in 

Knowing the density of states, we can calculate the partition function and its derivatives for 
any real or complex value of j3 . In particular, it can be used to determine the Fisher's zeros of the 
partition function ^]. Locating these zeros in the complex jS plane and their volume dependence 
is important to understand the large order behavior of the weak coupling expansion [Q, ^ |T^, 1 1 ] 
at zero temperature and the nature of the finite temperature transition [12|. Related questions have 
also been discussed in a poster presented at the same conference [13]. 



2. The density of states 

In the following, we focus on a SU (2) gauge theory with Wilson's action on a lattice and 
periodic boundary conditions. We denote the number of plaquettes = 6 x L^. The partition 
function Z(j8) is the Laplace transform of n{S), the density of states: 

Z{I5)= / dSn{S)e-^^ , (2.1) 
Jo 

with 

n{S)=Y[[ dUi5{S-Y,{\-{\/N)ReTr{Up))) (2.2) 

We can interpret ln(«(5')) as a "color entropy" (extensive). For cubic lattices with an even number 
of sites in each direction and a gauge group that contains — 1, it is possible to change pReTrUp into 
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—pReTrUp by a change of variables Ui ^ —Ui on a set of links such that for any plaquette, exactly 
one link of the set belongs to that plaquette [|l^]. This implies 

Z(-/3)=e2'5'^^Z(j8) (2.3) 

and consequently 

n{2,yVp - S) = n{S) (2.4) 



Thanks to this symmetry, we only need to know n{S) for < 5 < ^ . Note that < S >= J/'p 
means < TrlJp >=0. 
We define 

f{x,.^p) = ln{n{x.A'p,.^p))/^p . (2.5) 



The symmetry ( [2.4[ ) implies that 

/(x,^p)=/(2-x,^p) (2.6) 
The existence of the infinite volume limit requires that 

lim,j^^^„f{x,J^p) = f{x) , (2.7) 



with f{x) volume independent. In the same limit, the integral ( [2. 1] ) can be evaluated by the saddle 
point method. The maximization of the integrand requires 

fix) = p . (2.8) 



3. The one plaquette case 

In the case of the one plaquette model, the density of state simply follows from the explicit 
form of the Haar measure: 

nip,.{S) = -VSi2-S) (3.1) 
K 

At leading order, the large j8 behavior of the partition function is determined by the behavior of 
n{S) near S = 0. The fact that n{S) oc y/s for small S implies Z oc jS^^/^ for large j3. The l/j3 
corrections can be calculated by expanding the remaining factor \/2 — S in powers of S. One then 
sees that a series with finite radius of convergence becomes an asymptotic series if we integrate 
over S from to oo (instead of to 2). In addition, the large order behavior of the asymptotic series 
is determined by the non-analyticity of n\pi {S) at the maximal value of S (2 in this case). 

These properties are in agreement with the general idea that the large order behavior of the 
weak coupling expansion is determined by the behavior at small negative coupling [|l^, |T^. In the 
present case, small negative means that j3 is very negative. In this limit, the largest possible 
values of S dominate the integral (in agreement with what we explained above). It would be inter- 
esting to understand if this property persists on Lf' lattices. Unfortunately, numerical values of the 
weak coupling expansion of the plaquette are not available for SU (2) and we will have to rely on a 
model proposed in [^]. 
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4. Approximate forms of n (5) 

Numerical calculations of n{S) can be obtained by patching plaquette distributions multiplied 
by the inverse Boltzmann weight at various values of j8. In we presented numerical data for 
lattices with L = 4, 6 and 8. For these values of L, finite volume effects are not too large 
and plaquette distributions are broad enough to allow a reasonably smooth patching. The volume 
dependence is resolvable only for small values of 5 where a behavior ln{S)/'V is observed for 
ln{n{S)). The coefficient of the singularity was calculated to be (3/4) — {5/\2)L^^ in reasonably 
good agreement with the numerical data. 

The numerical results for f{x) were compared with expansions that can be obtained from the 
strong and weak coupling expansions of the average plaquette. Intermediate orders in these ex- 
pansions show a good overlap for values of S that correspond to the crossover (see Fig. |l|). The 
convergence of the new series can be related empirically to those of the series for the average pla- 
quette. The general picture that was obtained by trying with known series is that the converted 
series inherits the asymptotic behavior of the original series. The conversion of the series is per- 
formed using the saddle point equation. For the strong coupling, we expand about x=\ (we remind 
that X = S/^p, see section 2). Graphs of the accuracy of the expansion at successive orders, show 
a crossing characteristic of a finite radius of convergence near x = 0.5. This is consistent with a 
crossing near j3 ~ 2 for the plaquette (for j8 = 2, the average plaquette is about 0.47). For the weak 
coupling, we expand about x = 0. Accuracy graphs show consistent improvement as the order in- 
creases (with possible saturation) when x < 0.4 for f{x) and j8 > 3 for the plaquette. However it 
should be kept in mind that the large order of the series for the plaquette has been modeled rather 
than calculated explictly. For details and graphs see [^]. 

The weak coupling expansion determines the logarithmic singularities of ln{n{S)) at both 
boundaries. When these singularities are subtracted we obtain a bell-shaped function that can 
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Figure 1: Weak and strong coupling expansion of / at a few intermediate orders. 
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be approximated by Legendre or Chebyshev polynomials. Empirically, the determination of the 
expansion coefficients based on the discrete orthogonality of the Chebyshev polynomials (rather 
then numerical interpolation followed by numerical integration) seems the most stable. 



5. Calculation of Fisher's zeros 

One motivation for this work is to improve our ability to determine the zeros of the partition 
function in the complex j8 plane. For reference, it is useful to understand the limitation of the 
reweighting MC method. In order to estimate the errors in the location of the complex zeros, we 
considered the changes in the location of the zeros of the real and imaginary part due to statistical 
fluctuations. We considered 200 sets of 40,000 values of S picked at random out of the large sample 
of values (bootstraps) generated for j3 = 2.225. For each of the 200 sets, we calculated the zeros of 
the real part on a small grid with typical distance between neighbooring points of the order of 10^^. 
Using this procedure, 383895 zeros of the real part were found. We then studied the distribution of 
these zeros using a 200 by 200 grid in the j3 complex plane. The results are shown in Fig. 2. In 
this contour plot, the outer contours go through the bins that have 20 zeros, the first inner contours 
correspond to 60 zeros, the next to 100 zeros etc.. The circle of confidence [§] in the Gaussian 
approximation for 40,000 independent configurations as well as another estimate (red boxes in Fig. 



2) of the region of confidence discussed in [ ]10| ] are shown on this graph for reference. It is clear 
that as we get closer to the boundary of the region of confidence, the distributions get wider. 

It is easier to look at horizontal sections of this distribution. We then have simple histograms 
with 200 bins. The results are shown in Fig. 3 for /m j3 = 0.1, 0.115, 0.13 and 0.145. This allows 
us to observe the broadening of the four central peaks as Im j3 increases. For instance, the two 
most central peaks are quite narrow up to /m j3 = 0. 1, but their width becomes comparable to 
their separation when Im ^ > 0.13. One should bare in mind that such distributions should be 
understood together with the interlaced distributions for the imaginary part which follow similar 
patterns. It is clear that complex zeros found in regions where there are broad distributions are 
unreliable. The improvement in this situation obtained by using the j8 independent density of states 
presented above is discussed in a poster [p^]. 



6. Conclusions 

We have calculated numerically the density of states for SU (2) lattice gauge theory. The 
intermediate orders in weak and strong coupling agree well in an overlapping region of action 
values as shown in Fig. |l]. However, the large order behaviors of these expansions appear to 
be similar to the corresponding ones for the plaquette. Volume effects can be resolved for small 
actions values. Corrections to the saddle point estimate need to be developed more systematically. 
Aprroximation of a subtracted quantity by Chebyshev polynomials looks very promising. We also 
plan to use this method to study abelian gauge theories and the large N behavior of SU {N) gauge 
theories where interesting results based on the density of states have already been obtained [p7|]. 
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Figure 2: Distribution of zeros of the real part of the partition function in the complex j3 plane and regions 
of confidence described in the text. 
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Figure 3: Horizontal sections of the previous graph described in the text. 
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